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ABSTRACT 

We have observed the 2-60 keV flux from N6C4151 change by a 

factor of two on a timescale of 1.5 days. We do not detect any fluctuations 

in excess of a factor of three on timescales less than four hours. During 

a total observation of -v- 11 days there was no statistically significant 

changes in spectral shape. The spectrum can be fit by a power law with 

22 

photon index a 1*42 + 0.06 and column density 7.5 + 0.5 x 10 
2 

at/cm . A 2a residual to this fit implies fluorescent Fe line emission 
with E.W. 240 eV. Both synchrotron self-Compton and thermal Compton 
models are consistent with the X-ray data. 
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I. INTRODUCTION 

N6C 4*1 51 has been the most Intensively studied Seyfert galaxy 
in both the X-ray and optical frequency bands. Medium energy, 2-15 
keV, proportional counter observations (Ives, Sanford and Penston 1976) 
combined with higher energy, 10 - 200 keV, scintillator results (Baity 
et ai. 1975; Paclesas, Mushotzky and Pelling 1977) have shown that the 
X-ray spectrum can be best described by a power law with a low energy 
turnover due to photoelectric absorption. However, there has not been 
a single observation which covers the 2-60 keV band. The combination 
of 3 Ariel 5 observations (Barret^, 1977) has Indicated that the continuum 
remained constant In both Intensity and spectral Index while the column 
density varied by a factor of 5. Elvis (1976) observed an increase In 
the 2 - 10 keV flux of '^1.7 ever a time scale of 3 days. He Interpreted 
this as due to a variation In the continuum, but the lack of spectral 
data does not allow the unambiguous determination of Its origin in either 
the continuum or the column density. 

We present here the results of an 11 day observation, days 145- 
156 of 1977, over the 2-60 keV band from the Goddard Space Flight Center 
Cosmic X-ray experiment on OSO-8, This experiment and the methods of 
analysis have been described by Serlemitsos et 1976. 

II. RESULTS 

A. Intensity Variability 

The 2-60 keV intensity of NGC 4151 (Figure 1) changed in an ir- 
regular fashion during our observation with the most obvious feature being 

a factor of 2 increase from .014 + .001 cts/cm sec to .027 .001 cts 

2 

cm sec over a rise time of 1.5 days. This rise corresponds to a change 
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In the 2 - 6 keV flux of 4.3 x 10’^^ ergs/cm^sec to 8.2 x 10‘^^ ergs/ 

2 

an sec using our best fit spectral parameters. The flux then decreased 

p 

to an approximately constant level of .02 cts/cm sec over the next 2-3 
days. We note that these changes are not sharp but are gr. .ual and can 
be approximated by exponentials with a rise time of 2 days and a decay 
tiTO of 'V 7 days. These flux levels are similar to the range of 4.8 - 
12.7 X 10"^^ ergs/cm^sec seen by UHURU (Ulmer 1977). Thus the entire 
range of reported variability can occur on a time scale of less than 1 
week. In contrast to the spectral changes accompaning the intensity vari- 
ability seen by Ariel 5 (Barr et ^ 1977) there was no change in either 
spectral index or column density seen during our observation with an upper 
limit of La .15 and aN^ ''-lx 10^^. This is quite similar to the intensity 
changes in Cen-A reported by Mushotzky et (1978) which had a time 
scale of -v-3 days and a 50% amplitude. 

There has been a report by Tananbaum et al_. (1978) of variability 
on time scales of 10 minutes or less up to a factor 6 2 

average rate. We have examined our data on time scales from 10 sec to 
3 hours and have found no evidence for such variations in the source 
inconsistent with fluctuations in the background. We can set a 3a upper 
limit of no fluctuation greater than a factor 3 on any time scale in excess 
of 10 seconds during our 11 days of observation. During our observation 
the source had the same intensity range as that seen in the period when 
Tananbaum et al. observed the strong fluctuations. We conclude therefore 
that these short scale fluctuations, if real, are much less frequent than 
would be implied by the UHURU data. 



- 4 - 


B. Spectra 

At all times during our observation the spectrum (Figure 2) can 

be adequately fit by a pov^r law with photon Index a = 1.42 + .06 and 

with absorption due to cold material In the line of sight with cosmic 

abundance (Brown and Gould 1970; Fireman 1974) of = (7.5 + 0.5) x 10^^ 

2 2 

cm with a x per degree of freedom^ 1.35 (Fig, 2a). A thermal fit 
Is not weil determined but we require kT ^ 70 keV, corresponding to a 
temperature at which the knee Is higher than our highest observing energy. 
The observation of the source at this time out to 200 keV (Beall et al. 

1978) requires kT > 150 keV. 

There is an Indication of a 6.4 keV Fe fluorescent line at the 
2.5 level. The best fit model has an equivalent width of 240 eV 
(lo error). If real, this indicates (as does the Ariel 5 edge data, Barr 
et jn_) a Fe overabundance of '^2 relative to a solar abundance of 4 x 10"^, 
The OSO-8 edge data are also consistent with this. Because of the extreme 
weakness of this feature we cannot comment on any variability. 

III. DISCUSSION 

Our results appear to confirm the original somewhat weaker result 
of Elvis that N6C 4151 varies on time scales less than a week. The time- 
scales seen by us and by Elvis are similar and perhaps indicate the scale 
size of the source of a few light days. Our absorption column density 
differs from any of those reported by Barr et from Ariel 5 and therefore 
confirms that the absorption in N6C 4151 varies on a timescale of months. 
However all spectral data are consistent with a constant spectral index. 

We shall try to construct a self consistent model of the nuclear 
region of NGC 4151 which can, in a qualitative way, explain the X-ray 
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absorption variations, the X>ray continuum and the optical and IR spectra 
of tiie source. 

1) Lw energy absorption, optical reddening and IR spectrum. 

In principle the X-ray column density should be linearly related 
to the optical reddening (Sorenstein 1975). However In NGC 4151 the relation 
found by Gorenstein to apply to material In our galaxy seems not to apply. 

The lowest X-ray column density seen, 3 x 10*^ at/cm^, would Imply a red- 
dening Ay Z 19 mag. However the limit on the optical reddening is Ay Z 
,2 mag (Wu and Wee<tnan 1978). We can interpret this as caused by a low 
dust to gas ratio since the reddening is due to dust while the X-ray column 
is due to gas. The required dust to gas ratio is'v 1.6 x 10 by mass 
or about 1/200^^ that of '‘normal". However, this is quite similar to 
that seen in galactic HI I regions (Gillett et al_. 1975). 

It has been shown for several galaxies (M82, NGC253) that the 
bulk of the IR flux is due to cold dust. Given the above constraints 
on dust surrounding the nucleus of NGC 4151 it is of interest to see if 
the lOu IN nuciear component could be due to dust associated with the 
X-ray absorbing material. Using 7 x 10 and a gas to dust ratio 
of 6 X 10^ implies a dust column density of 2 x 10“® gm/cm^. If the dust 
is optically thin and is silicate material we can calculate the required 
amount of dust to account for the observed lOv; flux (Beall ^ iX. 1978). 

If we assume a dust temperature T Z 200°K which is similar to that of 
HII regions in our galaxy we require a mass of dust ^ 36 in a spherical 
cloud 30 pc (or .3" at 20 Mpc) in size. This physically oversimplified 
model would account for the lOw flux but could not account for the 3.4y 
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flux because of Its low tenperature. Since the 3.4 flux Is variable 
on a timescale of less than a year (Penston et §X. 1974) it must be 
due to another component. If it Is thermal it iraist be hotter and there- 
fore smaller. 

Another possible model to explain the high but low In l«5C 4151 
would require most of the optical emission to originate In a region 
outside of the place where the X-reo' absorption occurs. This would explain 
why the optical light Is not strongly reddened by the large amount of 
dust Implied by the X-ray column density If normal dust to gas ratios 
are assumed. In this case one would not expect and A^ to be related. 

2) Implications of X-ray Absorption on Optical Observations 

In this section we show how the variability in X-ray absorption 
determined from Ariel 5 and 0^-8 Is consistent with optical observations. 

Models accounting for the optical emission line spectrun in 
N6C 4151 require at least two regions, a broad line region (BLR) producing 
the broad Balmer lines and a narrow line region (NLR) producing the narrow 
core of the Balmer lines and the forbidden lines. Since the X-ray emission 
presumably shines through both of these regions their properties should 
be reflected in the X-ray absorption and variations therein. 

The NLR is characterized by a size of ''-50 pc- (Ulrich 1973) and 

a density 2 x 10^ cm”^ (Boksenberg et al[, 1975) and a filling factor 

-3 -2 

(the volumne of space taken up by gas) e''<10 - 10 . This would give 

21 22 

X-ray column densities in the range 3 x 10 - 3 x 10 . Thus it seems 

probable that NGC 4151 would always show X-ray absorption corresponding 
to the material in the NLR. However there is not enough material in the 
NLR to account for the total X-ray absorption and, since the lines in 
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the NLR are constant. It Is difficult to see how this region could account 
for variation In the X-ray absorption. We shall, therefore, assume (Sargent 
1973, Barr et al* 1978) that most of the X-ray absorption Is due to 
material In the BLR. 

Typical models of the BLR have densities, n^ ^ 10^ cm"^, sizes 
1 X 10^^ cm (Anderson 1973) and filling factors of 10’^ - 10’®. Thus 

pi pA o 

one expects X-ray column densities In the 10 - 10 at cm range. 

The fact that the X-ray absorption varies on timescales of years or less 

Implies that we see only a few clouds in the line of sight at any one 

22 

time and that each cloud has -v- 5 x 10 . We can set limits on the 
ratio of filament size, r, to tranverse velocity, Vj, of these filaments 
such that r/Vj < 1 x 10^ sec, the variability time scale seen by Ariel 5. We 
know that these filaments must cover the X-ray source and thus must have 
a size r< 3 light days (2.5 x 10^® cm). Thus Vy ^ 2.5 x 10® cm/sec 
('^.Olc) which Is on the order of the Doppler width of the broad lines, 
a similar such model was suggested by Ives et 1976. 

If we assun« (Osterbrock 1978) that such velocities are circular 
the mass of material Inward of 2 light weeks (4 x 10^® cm) necessary is 

p 7 

Me = r Vj/gZ 2 X 10' This Is indirect evidence for large mass con- 
centrations in the core of Seyfert galaxies. If the BLR Is larger the 

3 

mass Is larger by r . 

3) The Mechanism of X-Ray Emission 

We shall assume for this discussion that the X-ray emission 
Is due to a non-thermal process which Is related. In some fashion, to 
the non-thermal optical and UV emission. We shall discuss the synchrotron 
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self-Con^ton (SSC) inechanistn (Jones* O'Dell and Stein 1974) and a Comptonized 
accretion disk model due to Katz (1976). 

a) In a SSC model m shall attribute the observed rise time 
of the X-ray flux to the filling of the emission region by an Injection 
of relativistic particles. We shall attribute the decline of the X-ray 

©nlsslon to a cont)1 nation of synchrotron losses* Compton losses and adiabatic 
expansion. 

In a SSC model the X-ray emission Is due to the Compton scat- 
tering of the photons produced by synchrotron emission off the same electrons 
that produced them. In this model the observables necessary to define 
a physical model are the ratio of Compton to synchrotron flux the 
slope of the synchrotron component a* and 2 parameters of the synchrotron 

emission alone* the turnover frequency due to synchrotron self absorption 

w%and the flux ai this turnover F , Given these observables we are then 
^ ''n 

able* using a series of relations developed by Jones, O'Dell and Stein 

(1974) to derive physical parameters such as the size, e^, the magnetic 

field B and the loss times t due to synchrotron losses and t to Compton losses. 

D V 

Because in NGC 4151 we have no direct evidence of a compact high 
frequency component the problem is not well determined. It is however of 
interest to see if a model which is consistent with the IR data and the micro- 
wave upper limits can be constructed. For this self consistent model we have 
B ^ 1.5 gauss, % 1.5 x 10^^ Hz, Fy^ 'v 4 Jy, = 1.1 x lO”^^ rads 2 light 

days), E^^ 2 x 10 ^ a synchrotron lifetime t^ 30 days and a Compton lifetime 

t ''' 1 day. Such a model has a break in the synchrotron spectrum due 
to energy losses at Vg 8 x 10^^ Hz (% 40u) and would thus have a steeper 
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slope in the IR than in the millimeter. The X-ray spectrwn should deviate 

2 2 

from a power law at yn vg /v^^, whore Vp is the relativistic factor 

of W»e electrons emitting at (Burbidge» Jones and O'Dell 1974) » in 
this nK>del ^ 200 keV, The first order Compton scattered flux would 
be the dominant continuum flux at E ^ 200 eV, The contribution of this 
component to the near IR, optical and UV flux is difricult to calculate 
because the short lifetimes of the electrons indicate that repeated injection 
or reacceleration is important. However in the high frequency limit the 
srpe of the power law should approach a' ^4/30 + 1. This slope of 
^1.45 seems to be near the range of values a''^l.l - 1.6 seen by Puschell 
and Stein (1978). We feel that the UV data of Wu and Weedman (1978) is 
consistent with a slope «' ^<1.4 as are the U,B,V data of Beall (1978). 

We can crudely estimate the optical variability in this model by 
asswning that the changes in flux are due only to changes in the number 
of relativistic particles. In this case the optical flux would vary by 
"vl.U for a factor 2 change in X-ray flux. Thus on a short timescale a 
factor 2 change in X-ray flux would be accompanied by Am^ -v . 15 , We also 
note that the Compton UV component adds an additional source of ionizati-'n 
and heating which should be included in models of the BLR (Mac Alpine 
1974). 

b) It has been suggested (Katz 1976; Schnopper et 1977) that 
a model in which an unexplained source of soft photons is Compton scattered 
by collisions with electrons in a hot cloud may explain the non-thermal 
emission from QSO's and Seyferts. A similar model has been advanced to 
explain X-ra/ emission from Cyg X-1 (Shapiro, Lightman and Eardly 1976). 
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Fitting the observed ti tries cal e cf variability and spectral Index to Katz's 

model we would derive a cloud size R 2 x 10^^ cm {,77 4,d.) an electron 

9 -3 

density n^ 2 x 10 cm and a mass necessary to constrain the cloud 
of M 1 X 10® This model gives the proper X-ray luminosity and agrees 
with the fact that Intensity variations are not accompanied by spectral 
variations. In this model the optical variability should be similar to 
the X-ray variability and thus factor 2 changes in X-ray flux should be 
Am^ .3 mag. 

IV. CONCLUSIONS 

The observational data show that N6C 4151 continuum va»*1es by at 
least a factor 2 on timescales less than 2 days with no change in spectral 
properties during this variation. The absorption variation indicated 
from a combination of Ariel 5 and OSO-8 data is consistent with the variable 
X-ray absorption due to moving clouds In the broad line emission region. 

The X-ray emission is consistent with both the synchrotron self Compton 
mechanism and thermal Compton models. Our data Imply that the X-ray emis- 
sion from N6C 4151 ccnnes from a physical region of the same order as 
the non-thermal compact cores of radio galaxies (e.g, CsnA). 


COUNTS /CM -SEC 














FIKIRE CAPTIONS 


Figure 1 


Figure 2a - 


Figure 2b - 


Light curve of N6C4151 in May- June 1977. The data are 
binned in 4 hour intervals. 

The pulse height spectrum of N6C4151. The solid line 
is the best fitting power law with absorption convolved 
through the detector response. Note the excess flux 
over the model in the channels near 6.4 keV. The insert 
shows the 90% confidence contour for this fit. 

The best fitting photon spectrimi for a fit without an 
Fe line. This differs fran 2a in that a model is necessary 
to invert the PHA spectrum. 


REFERENCES 


Anderson, K. S. 1970, J. 162, 743. 

Barr, P. , White, N.E., Sanford, P. W. , and Ives, J. C. 1977, M.N.R.A.S. 
181 . 43. 

Baity, W. A., Jones, T. W. , Wheaton, W. A., and Peterson, L. E. 1975, 
( Letters) 199 , L5. 

Beall, J. H., Rose, W. K., Dennis, B. R., Crannell, C. J., Dolan, J. F., 
Frost, K. J., and Orwig, L. E. 1978, B.A.A.S. 9. 558. 

Boksenberg, A., Shortrldge, K. , Allen, D. A., Fosbury, R.A.E., Pens ton, 

M. V.. and Savage, A. 1975, M.N.R.A.S. 173 . 381. 

Brown, R. L. and Gould, R. J. 1970, Phys . Rev. Dl , 2252. 

Burbldge, 6. R., Jones, T. W. and O'Dell, S. L. 1974, J. 193 . 43. 

Elvis, M. 1976, M.N.R.A.S. 177, 7. 

Fireman, E. L. 1974, J. 187 , 57. 

Glllett, F. C., Forrest, W. J., Merrill, K. M., Capps, R. W., and Soifer, 
B. T. 1975, J. 2m, 609. 

Gorenstein, P. 1975, A£. 198, 95, 

Ives, J. C., Sanford, P. W. and Penston 1976, J^. ( Letters ) 207 . L159. 

Jones, T, W, , O'Dell, S. L. and Stein, W, A. 1974, Ag^, ;J. 188 , 353. 

Katz, J. 1972, A£. J. 206, 910. 

MacAlpine, G. 1974, J. 193, 37. 

Mushotzky, R. F., Serlemitsos, P. J., Becker, R. H. , Boldt, E. A., and 
Holt, S. S. 1978, J. 220, 790. 

Osterbrock, D. E. 1978, preprint. 

Paciesas, W. S,, Mushotzky, R. F., and Pelllng, R. M. 1977, M.N.R.A.S, 



Penston» M. V., Penston, M. J., Selmes, R. A., Becklin, E. E., and 
Neugebauer. G. 1974, M.N.R.A.S. 169 , 357. 

Puschell, J.» and Stein, W. A, 1978, private communication. 

Sargent, W, L. 1973, in X-Ray and Gamma Ray Astroncmiy , Bradt, H. and 
Glacconi, R., eds. 

Schnopper, H. W., Epstein, A,, Delvaille, J. P., Tucker, W., Doxsey, R., 
and Oernigan, G. 1977, ( Letters) 215 , L7. 

Serlemitsos, P. J., Becker, R. H., Boldt, E. A., Holt, S. S., Pravdo, 

S. H., Rothschild, R. E., and Swank, J. H. 1976, Proceedings on a 
Symposium on X-Ray Binaries, NASA SP-389, p. 67. 

Shapiro, S. L., Lightman, A. P., Eardley, D. M, 1976, J 204, 187. 

Tananbaum, H., Peters, G., Forman, W. , Glacconi, R. and Jones, C. 1978, 
submitted to J. 

Ulmer, M. 1977, J. ( Letters ) 218 , 1. 

Ulrich, M. H. 1973, J. 1^, 51. 

Wu, C. C. and Weedman, 0. 1978, preprint. 



BIBLIOGRAPHIC DATA SHEET 


1 . Report No. 
TM 79551 

2. Government Aw»$sion No. 

3. Recipient's Catalog No. 

4. Title end Subtitle 

X-Ray Observations of a Flare in NGC4151 
from OSO-8 

5. Report Date 
June 1978 

6. Performing Organization Code 
661 

7. Author(s)R.F. Mushotzky(NAS/NRC) . S.S. Holt, and 
P.J. Serlemitsos 

8. Performing Organization Report No. 

9. Performing Organization Name and Address 

Code 661 

Cosmic Radiations Branch 
Laboratory for High Energy Astrophysics 
NASA/Goddard Space Flight Center 
Greenhellf MarylancL 2Q771 

10. Work Unit No. 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

12. Sponsoring Agency Name and Address 

14. Sponsoring Agency Code 

15. Supplementary Notes 

To be published in the Astrophysical Journal (Letters) 

16. Abstract 

We have observed the 2-60 keV flux from N6C4151 change by a factor of two 
on a timescale of 1.5 days. We do not detect any fluctuations in excess of 
a factor of three on timescales less than four hours. During a total obser- 
vation of 'V 11 days there was no statistically significant changes in 
spectral shape. The spectrum can be fit by a power law with photon index 
a 1.42 + 0.06 and column density Njj 7.5 0.5 x lO^^ at/cm2, A 2a 

residual to this fit implies fluorescent Fe line emission with E.W. -v 240 eV. 
Both synchrotron self-Compton and thermal Compton models are consistent with 
the X-ray data. 

17. Key Words (Selected by Author(s)) 

18. Distribution Statement 

19. Secdrity Classif, (of this report) 

u 

20, Security Classit. (of this page) 21. No. ut Patjes 72. Pure* 

U 15 


•For sale hy National Technical Informatics) S(Mvice, Spnnyjfield, Vngnija 22151 


osrc 44 

















